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Abstract 

We discuss the system-specific optimization of 
long-range separated density functional theory 
(DFT) for the prediction of electronic properties 
relevant for a photocatalytic cycle based on an 
Ir(III) photosensitizer (IrPS). Special attention is 
paid to the charge-transfer properties, which are of 
key importance for the photoexcitation dynamics, 
but and cannot be correctly described by means of 
conventional DFT. The optimization of the range- 
separation parameter using the ASCF method is 
discussed for IrPS including its derivatives and 
complexes with electron donors and acceptors 
used in photocatalytic hydrogen production. Par¬ 
ticular attention is paid to the problems arising for 
a description of medium effects by means of a po¬ 
larizable continuum model. 

1 Introduction 

Transition metal (TM) organometallic complexes 
have found a wide range of catalytic, medical, and 
biological applications. In view of the growing de¬ 
mand in energy, one of the most perspective appli¬ 
cations of TM complexes is their role in systems 
for conversion and storage of solar light energy 
into chemical form. Here, a variety of homoge¬ 
neous and heterogeneous schemes have been sug- 
gested.^i^ TM complexes are attractive because 
of their notable spin-orbit coupling. This facili¬ 
tates the absorption of sun light, whose energy can 
be stored in charge-separated triplet states. These 
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long-living triplet states are then available for fur¬ 
ther reactions. Besides light-harvesting properties, 
TM complexes are also used as direct catalysts 
for water splitting; for recent reviews on homoge¬ 
neous photocatalysis see, e.g. refs SHI, 

To develop new efficient and stable photocat¬ 
alytic systems, understanding of their primary 
photoreaction steps and especially excited state 
properties is required.!^ At the moment, most com¬ 
putationally demanding studies of sizeable TM 
complexes applied in catalysis are performed with 
DFT in combination with the B3LYP functional 
and its extension in the time-domain in the lin¬ 
ear response formulation (TDDFT).l2Hi^ Besides 
efficiency it is the absence of system-dependent 
parameters that need to be determined first (such 
as the active space in multi-reference methods), 
which makes (TD)DFT attractive. The power of 
the DFT method to reproduce different electronic 
ground state properties even for rather large sys¬ 
tems is well documented.® Concerning excited 
state properties of TM complexes, however, there 
appears to be no unequivocal opinion; for a review 
see ref®. For an Iridium (III) heteroleptic com¬ 
plex, for instance, we have shown that the excited 
state energies are strongly dependent on the em¬ 
ployed exchange-correlation functional.® Never¬ 
theless, the TDDFT approach remains to be the 
most attractive one. 

Many of the practical problems are related to 
the er roneo us description of charge transfer (CT) 
states.®® To overcome this drawback, different 
schem es hav e been proposed, such as, e.g., scaled 
hy brids® ® and range-separated hybrid function¬ 
als.®® 
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Benchmarking DFT results for TM complexes 
in solution is hampered by the facts that there 
are no high-resolution experimental data avail¬ 
able and only a limited number of theoretical ap¬ 
proaches might be applicable to generate a pre¬ 
sumably more accurate reference. At the mo¬ 
ment, multi-reference perturbation theory is the 
best choice for reference calculations of excited 
state properties. However, in multi-reference ap¬ 
proaches the problem of a sufficiently large and 
balanced active space has to be addressed for ev¬ 
ery system.^ 

In this paper, we focus on the tuning of the 
range separation parameter for the long-range cor¬ 
rected LC-BLYP functional as applied for a par- 
ticular photocatalytic system developed by Beller 
et al.l23 it consists of a heteroleptic Ir(III) photo¬ 
sensitizer, triethylamine (TEA) as a sacrificial re- 
ductant and a series of iron carbonyls (Fe-cat) as 
water reduction catalyst (Figure 1). In addition, 
hybrid systems of IrPS and silver nanoparticles 
are considered.l^^^ Finally, the scope is broad¬ 
ened by including results on modified IrPSs and 
copper based PS which could replace noble-metal 
containing PS, see refill and references therein. 
Although the choice of the system is arbitrary, its 
general properties as well as the basic reaction 
steps are common for TM photocatalysis. There¬ 
fore, the present study provides some general in¬ 
sight into the role of tuning the range separation 
parameter in these systems. 

In the following section, we start with a brief 
introduction in the theory of long-range separated 
functionals. We then present our results on tuning 
the range-separation functional for IrPS and joint 
IrPS-X systems, where X is a reactant responsible 
for a particular process, i.e. light absorption, elec¬ 
tron transfer and recovery of the ground state. The 
influence of the reference geometry, environmen¬ 
tal and basis set effects are discussed. In particular, 
the applicability of the ©-tuning together with the 
Polarizable Continuum Model (PCM) is critically 
analyzed. 


2 Optimally Timed Range- 
Separated Functionals 

DFT has its form al roo ts in the theorems of Ho- 
henberg and Kohnl^SEH and it can be considered as 
an in principle exact approach. In practice, it is 
hampered by the lack of the exact exchange cor¬ 
relation (XC) functional, Exc[p{^)]^ with p(r) be¬ 
ing the electron density. There seems to be no first 
principle route to Exc [p (i^)] and therefore it is usu¬ 
ally constructed on the basis of model systems or 
fitted to experimental reference data. Excited state 
calculations performed by TDDFT in the linear re¬ 
sponse approximation employ XC functionals ob¬ 
tained for the electronic ground state. 

All caveats and the sometimes spurious be¬ 
haviour of the DFT approach originate from the 
approximation of the unknown XC kernel. There 
are three main problems: First, in the r ^ oo limit 
the potential and hence the density itself possess 
an incorrect behaviour. In principle, at large dis¬ 
tances the potential should be dominated by the 
exchange term decaying as — 1 / r, whil e the cor¬ 
relation term decays as ~ _i^^4 |32l33] However, 
for the local density approximation and to a lesser 
extent for the generalized gradient approximation 
(GGA), the exchange term decays exponentially .1^ 
Second, the approximate nature of the exchange 
term does not cancel the Coulomb interaction of 
the electron with itself at large r as it is the case 
for the Hartree-Fock (HF) method (self-interaction 
error). Third, a well-known deficiency of ap¬ 
proximate Kohn-Sham approach is that the funda¬ 
mental gap, defined as a difference between ion¬ 
ization potential (IP) and electron affinity (FA), 
differs notably from the orbital energy difference 
£homo “ £lumo- ^ This can be explained by the 
finite jump of the Kohn-Sham correlation potential 
for a statistic ensemble with variable number of 
electrons while passing through integer (N) num¬ 
ber of electrons, the so-called derivative disconti¬ 
nuity. According to Koopmans’ theorem for 
the DFT case,l22I21ESlSl^he HOMO corresponds to 
the IP, but the LUMO is generally more strongly 
bound than in HF theory and cannot be related to 
the FA. 

Considering TDDFT calculations of electronic 
excitation energies, the error correlates with the 
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Figure 1: Scheme of photocatalytical water splitting introduced in ref 1^3 and structural formulas of pho¬ 
tosensitizers and water reduction catalysts studied in this work. 


overlap of the donor and acceptor orbit als, be ing 
the largest for CT and Rydberg states. 03^ For 
vanishing overlap between these orbitals (e.g., due 
to long-range CT), the excitation energy is reduced 
to the orbital energy difference, being a poor esti¬ 
mate in case of DFT, which is in contrast to HF 
theory. 

The remedy for the erroneous exchange poten¬ 
tial could be the substitution of the approximate 
density-dependent exchange energy with the exact 
orbital-depe ndent on el^ within generalized Kohn- 
Sham theory.l^^EIMl'pj^g exact exchange can be in¬ 
cluded in a fixed manner like in hybrid functionals 
(e.g. B3LYP) or weighted with a function depend¬ 
ing on the inter-electron distance ri 2 . The latter 
case is implemented in range-separated function¬ 
als via splitting the Coulomb operator into local 


and non-local parts 

1 ^ l-r(cQri2) ^ r(mri2) 

^12 ^12 ^12 

where r(u)ri 2 ) is a smooth range-separation func¬ 
tion, which damps the exchange contribution from 
the density functional and complements it with 
exact exchange. Examples are the Yukawa ker- 
nelM c-“7r or the error functionl^®^ kernel 
erf(©r)/r. Such an approach eliminates the spuri¬ 
ous behaviour because the exact exchange has the 
correct asymptotic character and cancels the self¬ 
interaction exactly.l^^^^ 

Certain standard values for the range-separ ation 
parameter co have been established in refs 1^213] 
and they are used like universal c onstan ts in pop¬ 
ular quantum chemical programs .^*^1 The value 
0.33bohr“^ was determined by a least square fit 
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to empirical data for first- to third-row atoms^^ 
and later refined to 0.47bohr“^ for larger sets 
of small molecule s.l^ (for other test sets, see 
also refs.l^2l50l56|57l) However, it is clear that for 
molecules, which have not been part of the training 
sets, a more accurate description will be provided 
by choosing a system-dependent to. An optimal co 
for a particular system can be determined by en¬ 
suring that the energy of the HOMO orbital equals 
the negative of the IP, a relation that would be ful¬ 
filled for the exact functional.l^^ 

However, in the following we will u se alternative 
the so-called ASCF method,l^2lS9]60| t^g jp 

and EA are calculated as the differences between 
ground state (gs) energies of systems with N and 
A ± 1 electrons, i.e. 

IP“(A) = Eg®(A- 1) -Eg®(A), (2) 

IP®(A+ 1) = EA®(A) = Eg(A) -Eg®(A+ 1). 

(3) 

This yields the separate tuning conditions 

^(®) = |eHOMoW + IP“WI. w 


A(®) = I<!homo('V+1)+EA‘»(W)|. (5) 

In order to obtain a proper description of the fun¬ 
damental gap, the functions Jq(Q)) and Ji(co) for 
IP and EA should be minimized simultaneously. 
For the present one-parameter formulation this re¬ 
quires to minimize the general function 

/(m) =/o(©)+/i(©). (6) 

Note that in general /(to) can be non-zero even for 
exact functionals, because it is defined for systems 
with different numbers of electrons. 

Although minimizing /(co) was applied exten¬ 
sively for tuning the range-separation parameter, 
in ref it was stressed that the m-values are very 
close to minima of Jq(co), eq 4, or 7j(®), eq 5, 
depending on the particular form of the dependen¬ 
cies of £homo energy on ®. To avoid 

such a biased behaviour, it was suggested to apply 
a least-square approach to minimize the resulting 


error according 

r{0i) = yJjl{0i)EJl{(0). (7) 

Alternatively, for long-range CT systems where 
donor and acceptor units can be clearly distin¬ 
guished, ground states of neutral donor and nega¬ 
tively ionize d acce ptor can be considered in a tun¬ 
ing approachl^SES based on Mulliken’s rule.^^ 

It is important to stress that the optimal (O is 
density and thus system dependent. Thereby, 
the dependence of (O on the density for range- 
separated functionals is much more pronounced 
than that of the scaling factor for hybrid function- 
als.l23 Since the range- separa ted functionals elim¬ 
inate self-interactionBS^i^ and have the correct 
asymptotic form of non-local exact exchange, the 
correct Coulomb-like behaviour of the energy is 
assured and thus CT and Rydberg excitation ener¬ 
gies are improved. Additionally, since the tuned 
range-separation parameter improves the funda¬ 
mental gap for the A electron system, one can ex¬ 
pect also a better description of electronic excita¬ 
tion energies in TDDFT, where the leading term 
reads as Chomo — £lumo- 

In the past it was shown in numerous publi¬ 
cations the the optimally tuned functionals in¬ 
deed improve the description of different molec¬ 
ular properties such as IPs, fundamental and op¬ 
tical gaps, and CT and Rydberg transition en- 

ergics. ™35i40i42i50i52i53i59t l6Tfell^ accuracy of 

optimally tuned range-separated functionals was 
critically tested in ref^^ for basis set variation 
as well as prediction of relative energies of spin 
states, binding energies, and the form of potential 
energy surfaces. 

Finally, we note that a generalized form of LC- 
partitioning of Coulomb interaction proposed by 
Yanai for CAM-B3LYF® is often applied for 
non -empi rical optimal tuning as well (see, e.g. 
refs F6|79H 83l^ Additionally, the simultaneous two- 
parameter optimization of short- and long-range 
separation parameters was proposed in order to 
minimize the delocalization error. Successfully 
applications provide d an imp roved description of 
molecular properties. 
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Figure 2: Overview of all optimized structures of those molecules that were used for ©-tuning. 


3 Computational Details 

The tuning of the range-separation parameter (O 
according to Eqs 6 and 7 h as been done for the 
LC-BLYP functional.l ‘^^ l ^^l^‘^ 1 It includes the Becke 
exchangel^and a correlation pa rt of G GA-type es¬ 
tablished by Lee, Yang and Parr.l2^*^This scheme 
uses eq 1 with an error function kernel as sug¬ 
gested by Hirao et. ak.l^ For the case of the par¬ 
ent IrPS compound, the results with a tuned (O 
are compared to those obtained with the standard 
system-independent values of 0.33 bohr“^ and 
0.47bohr“^ as well as with BLYP and B3LYF® 
functionals. 

All calculations have been performed using the 
LANL2DZ ECP basis set for Ir, Fe, Cu, and Ag 
and the 6-31G(d) basis set for all other atoms. In 
ref El it was pointed to a strong dependence of the 
optimal (0 on the inclusion of diffuse function in 
case of simple atomic and diatomic systems. To 
test the influence of diffuse basis functions for the 
present example of IrPS, calculations have been 
performed with 6-31-i-G(d) and 6-31-i-i-G(d) basis 
sets. The inclusion of one or two diffuse function 
has led to decrease of (O by 0.01 bohr“^ Since 
this effect appears to be of minor importance, be¬ 
low only results for the 6-31G(d) basis set will be 
presented. 


The impact of reference geometry used for tun¬ 
ing of CO was studied on the example of IrPS. 
The tuning for the slightly different geometries, 
which have been obtained within LC-BLYP with 
the two standard (0.33 and 0.47) and the optimized 
(0.18) CO, led to the same value of 0.18bohr“^. In 
the following, the structures optimized with LC- 
BLYP with optimal CO have been utilised for tun¬ 
ing (see Figure 2). For the complexes IrPS—X 
(X=TEA or Fe-cat) the geometries of the consti- 
tuing parts have been first optimized separately 
and then placed at fixed positions without further 
optimization (for further details see ref^^. For 
systems containing silver clusters, the geometries 
were optimized in refsl2^^^. The standard TDDFT 
formalism was applied for excited state calcula¬ 
tions. Solvent effects have been included within 
the PCM approach.l^ All calculations were done 
with the Gaussian09 set of programs.l^ 

4 Results and Discussion 

4.1 Optimization of the Range- 
Separation Parameter 

The optimal values of the range-separation param¬ 
eter CO for molecules and complexes related to the 
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Figure 3: Left panel: The functions defined in eqs 4-7 used for optimization of the range-separation 
parameter co for IrPS—Ag 2 . Right panel: Selected examples of J{(0) for IrPS in different environments. 


studied photocatalytic system for hydrogen gener¬ 
ation (Figure 1) as well relevant literature data for 
similar substances are collected in Table 1. In the 
following, we will discuss particular aspects of the 
optimization for the chosen target systems. 

Exemplary, in the left panel of Figure 3 details 
of the tuning procedure for the IrPS—Ag 2 case are 
presented. Since Jq{(o) and Jiico) describe sys¬ 
tems with different numbers of electrons, in gen¬ 
eral, their minima do not coincide. The shown 
example is the only case among the studied sys¬ 
tems, where the minimum of J{(0) is not well- 
defined due to an almost flat region from 0.18 — 
0.25bohr“\ where Jq{o)) and Ji{(o) upon sum¬ 
mation compensate each other. If the least-square 
function J*{(o) is applied instead, a minimum at 
0.19bohr“^ can be easily located. For all other test 
systems collected in Table 1, the minima of J and 
J* were very close to each other; J* led in some 
cases to a decrease of to by 0.01 — 0.02bohr“^ 
Since the use of J* was not crucial in our case, we 
focused on the results obtained by applying /, as it 
was the main approach of previous investigations. 

In the right panel of Figure 3 selected exam¬ 
ples of /(to) for IrPS in different environments are 
presented. The curve for IrPS(vac) represents a 
case with regular behaviour, i.e. there is a dis¬ 
tinct minimum. For IrPS—AgjQ the region af¬ 
ter 0.25bohr“^ corresponds to a changed order of 
HOMO and HOMO-1 orbitals for N + 1-electron 
system. For the isolated reduced IrPS in the dou¬ 


blet electronic state, the unpaired electron was 
previously shown to be localized on 7r(bpy)* or- 
bital.l^For the IrPS —Ag^Q, the internal oxidation- 
reduction occurs for (O > 0.25bohr“^: the un¬ 
paired electron moves from ;r(bpy)* to (7*(Ag). 

As IrPS forms no stable complexes, IrPS—X, 
neither with TEA nor with the iron-catalysts^, 
the relative position of constituents might be flex¬ 
ibly varied. Similar to the procedure in refl^ 
we have varied the position of X=TEA or Ee-cat 
on a sphere around IrPS. Thereby we found that 
choosing different locations does not have an im¬ 
pact on resulting CO. However, changing the dis¬ 
tance between IrPS and X led to some variation 
of optimal co, see Table 1. Interestingly, the op¬ 
timal CO value was found to increase in the range 
0.23 — 0.32bohr“^ with increasing the distance be¬ 
tween IrPS and TEA from 7 to 12 A. While the 
minimal value of Jq{co) remains at 0.18bohr“\ 
the minimum of (to) occurs at larger co. The us¬ 
age of J*{co) instead of /(to) provides a narrower 
range of 0 .23 — 0.28bohr“^. In passing we note 
that in refs it was shown that an optimization 
of CO for each geometry and subsystem can lead to 
size-inconsistency or unphysical potential energy 
curves in excited states. 

Eor comparison, the standard^^^l^ system- 
independent = 0.33bohr“^ and ®2 = 
0.47bohr“^ are shown in both panels of Fig¬ 
ure 3 by vertical lines. As can be seen from 
that figure and Table 1, the optimally tuned co 
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Table 1: Optimized range-separation parame¬ 
ters, ®, defined by the minumum of 7(®), eq 
6, for different molecules depicted in Figure 
1. For comparison, literature data on molecules 
relevant for photovoltaics are provided. 


Compound 

Optimal CD [bohr 

IrPS 

0.18 

IrPS (THF, PCM) 

0.05 

IrPS -1- 28 THF 

0.15 

[Ir(ppy) 2 (en)]+ 

0.18 

[Ir(pox) 2 (bpy)]+ 

0.18 

[Ir(pbox) 2 (bpy)]+ 

0.18 

[ii'(ppy)2(bpym)ir(ppy)2]^ 

0.14 

[Cu(phen*)(xant)]+ 

0.16 

Ag2 

0.42 

Agio 

0.22 

Ag20 

?«0.20‘‘ 

IrPS-Ag 2 

0.18-0.25^ 

IrPS—Ag^Q 

0.16 

IrPS—Ag 2 o 

0.16 

TEA 

0.31 

IrPS-TEA 

0.23 - 0.32 

Fe-cat(l) 

0.24 

Fe-cat(2) 

0.19 

Fe-cat(3) 

0.17 

IrPS-Fe-cat(l) 

0.18" 

IrPS-Fe-cat(2) 

0.17-0.18" 

IrPS-Fe-cat(3) 

0.16-0.17" 

r’ llUl 

'“60 

0.21 

phthalocyaninel^ 

0.16 

Si-nanocristals (5-15 A)l^ 

0.10-0.24 

coumarine dyes^^ 

0.17-0.21 

pentacene-C^I^ 

0.21 

nucleobases^^Sl 

0.27-0.31 


‘‘Because of convergence problems, only 
tentative a value is given; ^ using 7(ft)), eq 
6, a flat dependence is obtained, see text and 
Figure 3; 

different geometries are used, see text. 


values for compounds relevant in photocataly¬ 
sis are substantially lower than those deter mined 
for diatomics and small molecules.Since 
(0~^ reflects a characteristic distance for switch¬ 
ing between short- and long-range parts or, in 
other words, an effective electron screening (de- 
localization) length, previously optimal (O val¬ 
ues were found to decrease wit h increasing sys- 
tem size and conjugation length.l^ ^^l^^ l ^^l^ l ^^ l ^^ l 
However, in some cases the dependence was not 
monotonous and strongly varied for systems with 
different electronic structure.!^ In the present 
study, the Ir(III) complexes have various sizes, 
being smallest for [Ir(ppy) 2 (en)]''' and largest for 
[Ir(ppy) 2 (bpym)Ir(ppy) 2 ]'‘'. However, only a mi¬ 
nor dependence of (O on the size of ligands around 
the central Ir atom and the size of molecules (sil¬ 
ver clusters, iron carbonyls) bound to the PS was 
found (cf. Table 1). 

4.2 Effect of Solvation 

By requiring Koopmans’ theorem to be satisfied 
simultaneously for systems with N and +1 elec¬ 
trons where HOMO and LUMO are on differ¬ 
ent fragments, we automatically improve the de¬ 
scription of electron transfer. This implies that 
the long-range corrected functional with optimized 
range-separation parameter should have a predic¬ 
tive power to study also redox properties. To de¬ 
scribe redox (electron transfer) reactions in pho¬ 
tocatalysis, solvents effects need to be included. 
However, to the best of our knowledge, so far no 
m-tuning in the presence of a solvent was reported. 
The green curve with square symbols in Figure 3 
shows /(to) for IrPS, solvated via the PCM model 
(THF solvent). This yields a very small opti¬ 
mal CD = 0.05bohr“^ Similar computations have 
been done for TEA and Fe-cat(3), resulting in val¬ 
ues 0.10 and 0.05 bohr“\ respectively. Note that 
such a small value of CD implies that there is al¬ 
most no long-range correction from exact HF ex¬ 
change. Interestingly, the inclusion of 28 explicit 
THF molecules around IrPS, which corresponds to 
one solvation shell, led to CD = 0.15bohr“^ similar 
to vacuum case. 

According to Janak’s theoreml^ the redox en¬ 
ergy upon changing the number of electrons from 
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Figure 4: ASCF redox energies AE (dashed lines) and orbital eigenvalues — ^homo/lumo (solid lines) for 
different charge states (+2,+l, and 0) of IrPS computed with pure BLYP and LC-BLYP with (O optimiza¬ 
tion in vacuum and PCM solvent. The eigenvalue of the /3 spin-orbital was used for LUMO(-i-2) all other 
values correspond to the a spin-orbital. 


N—1 to N can be obtained in integral form as 


N 



N-l 


In the absence of self-interaction and delocaliza¬ 
tion errors, the orbital energies e stay constant 
between integer numbers of electrons, thus satis¬ 
fying the condition of 0) optimization, eq 2. In 
Figure 4 redox and HOMO/LUMO energies are 
shown for different functionals and optimization 
conditions. For the tuned to, AE = —^nouoi^) — 
“^LUMo(^ — 1) are given in central panel of Fig¬ 
ure 4 for IrPS in vacuum. When using the conven¬ 
tional BLYP functional this condition is not satis¬ 
fied (left panel). However, if we include implicit 
PCM solvation the jump in energies of HOMO(A^) 
and LUMO(A^ — 1) is mitigated by the response 
of the polarizable continuum. This leads to much 
better estimates of redox energies obtained with 
conventional functionals with PCM as compared 
with the vacuum.l^ In other words, in terms of 
the derivative discontinuity the PCM solvent ef¬ 
fectively leads to a partial mitigation of the ^xc 


deficiencies. 

That is why one needs to include much less exact 
exchange into the LC-BLYP functional (smaller 
(o) to satisfy Koopmans’ theorem and derivative 
discontinuity condition as is seen from Figure 3 
and Table 1. The orbital eigenvalues in the right 
panel of Figure 4 for the solvent case are almost 
constant, the deviations being probably due to a 
residual delocalization error. Importantly, the 0) 
value optimized for the vacuum does not fulfil the 
condition of constant orbital energies in solvent 
and vice versa. Further examples of Fe-cat(3) and 
TEA leading to the same conclusion can be found 
in the Supporting Information. 

4.3 Triplet Stability 

The nature and energetic position of the lowest 
triplet state of PSs need to be calculated very ac¬ 
curately since this state plays a key role in photo¬ 
processes as a dominant long-living and emitting 
excited state. However, the calculation of triplet 
states represents a challenge for the TDDFT ap¬ 
proach with long-range corrected hybrid function- 
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Figure 5: Energies of the lowest five singlet and 
triplet excited states and the lowest eigenvalue of 
the ground state stability matrix (A) in dependence 
on the range separation parameter co for IrPS in 
vacuum. 

als because of symmetr y brea king instabilities of 
the ground state solution .12512^ Typically, the larger 
amount of exact exchange leads to more pro¬ 
nounced problems with instabilities of the ground 
state wave function. This results in a divergence 
of TDDFT excitation energies and an imaginary 
energy of the lowest triplet state.l^ As shown in 
refl221, the non-empirical tuning of a> could avoid 
the instabilities through the stabilizing the ground 
state solution, but this issue should be carefully an¬ 
alyzed for the particular system under study. Since 
the instability problems are generally m ore p ro- 
nounced for the HF method than for DFT,I^2125| 
application of a smaller value for to should result 
in a more stable ground-state solution. In Figure 
5, the dependency of the smallest eigenvalue of 
triplet stability matrix!^ on the range-separation 
parameter, A(©), for the case of IrPS is given; the 
energies of lowest triplet and singlet energies are 
provided as well. 

It can be clearly seen that for the present sys¬ 
tem the maximum of stability corresponds 
to the optimal (O value. The A va lue is notably 
higher for the standard values © 12 - ^^^^This im¬ 
plies that LC-BLYP with properly tuned range- 
separation parameter ensures a more reliable es¬ 
timate of triplet state energies. The singlet state 
energies systematically rise with the increase of 


the exact exchange part, which is a co m mon be¬ 
haviour for TDDFT.I^For triplet states, this rising 
is compensated by the decrease of energies upon 
increase of instability resulting in maxima around 
0.25 - 0.40 bohr-^ 

4.4 Electronic Absorption Spectra 

A comparative analysis of excitation and differen¬ 
tial absorption spectra of IrPS in THF calculated 
by means of local (BLYP), hybrid (B3LYP), and 
long-range corrected LC-BLYP with standard and 
both optimal (for vacuum and THF) (O values is 
provided in Figure 6. In the TDDFT calculations, 
the 80 lowest singlet-singlet states have been in¬ 
cluded and assigned. A phenomenological broad¬ 
ening by a Gaussian lineshape with a width of 0.1 
eV was assumed. 

Overall, the features of the absorption spec¬ 
trum are well reproduced by B3LYP and LC- 
BLYP with optimal 0) = 0.18bohr“^. However, 
absorption and difference spectra could only be 
properly assigned simultaneously using LC-BLYP 
(© = 0.18bohr“^).l^ It should be noted, that be¬ 
sides the overall shape of spectra, which can be 
described by the conventional hybrid functional, 
the quality of description of CT states and band 
gap characteristics is only assured for the long- 
range functional with the properly tuned (O param¬ 
eter. The excessive amount of HF exchange (as in 
LC-BLYP with both standard 0) values) as well as 
complete absence of exact exchange (as in BLYP) 
led to shapes of absorption and difference spec¬ 
tra of IrPS deviating from experiments. Naturally, 
the results of LC-BLYP with optimal (O tuned for 
PCM are very close to those of pure BLYP be¬ 
cause of only small amount of exact exchange is 
included. Hence, LC-BLYP with © = 0.05bohr“^ 
optimized in a PCM solvent model does not agree 
with experiment, similar to the parent BLYP func¬ 
tional which is known for wrong prediction of en¬ 
ergies of CT states. 

The latter observation deserves some more com¬ 
ments. In view of the derivative discontinuity, © = 
0.05bohr“^ obtained for the PCM model gives an 
improved description as compared to the gas phase 
value (see Figure 4, right panel). In contrast, us¬ 
ing the gas phase © = 0.18bohr“^ together with a 
PCM model yields an “overscreening” (see Figure 
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Q 

O 




Figure 6: Left panel: Experimental (optical density, OD) and theoretical (TDDFT, oscillator strength) 
electronic excitation spectra. Right panel: Difference absorption spectra (reduced species - oxidized 
species) obtained theoretically and by transient absorption experiments of IrPS in THF-solution. Experi¬ 
mental spectra were measured by Loch brunner et al. and had been published together with the calculated 
BLYP and B3LYP spectra in refs^^^^^. All calculations included THF within the PCM. 


4, middle panel). At this point, one should note 
that under the influence of a polar solvent the en¬ 
ergies of CT and local states shift differently from 
their vacuum values because of substantially dif¬ 
ferent dipole moments. That is why, to correctly 
predict the electronic absorption spectrum in solu¬ 
tion, one needs to use to tuned in vacuum and to in¬ 
clude solvent effects via PCM. In other words, we 
observe that, when using the PCM model, for ob¬ 
taining a reasonable agreement between measured 
and calculated absorption spectra one has to sacri¬ 
fice the derivative discontinuity condition. 

5 Conclusions 

The present investigation of molecules and com¬ 
plexes relevant to photochemical water splitting 
points to the need for a system-specific optimiza¬ 
tion of the range-separation parameter within the 
LC-BLYP approach. It provides an improved 
band gap characteristics, increased ground state 
stability, and a correct asymptotic behaviour of 


CT states. The actual optimized value, (O = 
0.18bohr“\ was found to be nearly the same 
for various Ir and Cu-based PS as well as for 
IrPS complexes with silver clusters and iron car¬ 
bonyls. Using the CO optimization it was shown 
that the tuning for supermolecular compound sys¬ 
tems, where an electron relay between constituents 
is of interest, gives an improvement of the energies 
of HOMO and LUMO orbitals localized on differ¬ 
ent fragments. Thus, this approach is in principle 
suitable for accurate modelling of redox molecular 
properties. 

The optimization of the range-separation param¬ 
eter is usually performed for molecules in vac¬ 
uum. Here, for the first time, an optimization in 
the presence of a solvent described by the PCM 
model has been done. For the used ASCF approach 
this yielded some unexpected results, namely that 
the optimal to tends to zero. As a consequence 
there is almost no long-range correction and the 
quality, e.g. of electronic excitation spectra, cor¬ 
responds to that of the BLYP functional. Since 
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CT and long-range properties are of primary in¬ 
terest for photocatalysis, as a conclusion it is rec¬ 
ommended to use CO optimized in vacuum together 
with PCM. This seems to provide a better de¬ 
scription of CT processes and electronic absorp¬ 
tion spectra, even though the formal conditions im¬ 
posed by the derivative discontinuity are violated. 
Therefore, if redox reactions are studied, the com¬ 
bination ASCF/PCM is not recommended. It ei¬ 
ther leads to derivative discontinuity problems or 
to a wrong asymptotic behaviour of the optimized 
functional. Needless to say, that this procedure 
should not be applied without careful diagnostics, 
e.g., as provided by Figure 4. 
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Fig. SI. ASCF redox energies (AE) and orbital eigenvalues (— %omo/lumo) different charged states of 
IrPS, Fe-cat(3), and TEA computed with pure BLYP and LC-BLYP with different co in vacuum and PCM 
solvent (THF). Values of co are optimized for IrPS in vacuum (0.18bohr“^) as well as for each compound 
in vacuum and THE (given at the corresponding panels) 






































